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Many systems respond to slowly changing external conditions with crackling noise, created by 
avalanches or pulses of a broad range of sizes. Examples range from Barkhausen Noise in magnets to 
earthquakes. Here we discuss how the scaling behavior of the avalanche size and duration distribution 
and the power spectra of this noise depend on the rate O. at which the external conditions are 
changed. We derive an exponent inequality as a criteria for the relevance of adding a small driving 
rate f2 to the adiabatic model. We use the zero temperature nonequilibrium random field Ising 
model to test our results, which are expected to be applicable to a large class of systems with 
crackling noise. They also agree with recent experiments on Barkhausen noise in various materials. 

PACS numbers: 75.60.Ej, 05.40.-a, 64.60. Ht 



Jerky response to smoothly varying forces has been 
enjoying increasing attention for more than a decade. 
One reason for this is the surprising universality ob- 
served in the noise statistics of vastly different sys- 
tems. Earthquake faults, magnetic tapes, vortices in 
type II superconductors, charge density waves, mechani- 
cally failing materials, and many other systems respond 
to slowly changing external conditions with discrete, 
impulse-like events ("avalanches", or "pulses", leading to 
familiar "crackling noise" J|]) that span a huge range of 
sizes: earthquakes for example are known to range from 
tiny tremors to devastating magnitude 9 quakes on the 
Richter scale ||. Martensites crackle when subjected to 
slow temperature or stress changes [||, Q| . Similarly, when 
a magnetic tape is slowly magnetized by an external mag- 
netic field, the process happens in steps ("Barkhausen 
noise" (BN)), by avalanches of re-orienting magnetic do- 
mains, that range from microscopic to macroscopic in size 
H . The avalanche size and duration distributions are of- 
ten described by universal power laws over a broad range 
of sizes. The scaling behavior on long length scales is in- 
dependent of the microscopic details, and only dependent 
on a few basic properties, such as symmetries, dimen- 
sions, interaction range, and dynamics (and a few others). 
In many systems this universality has been explained the- 
oretically as being due to either proximity to an under- 
lying critical point [[|, or so called "self-organized criti- 
cality" (SOC) @. 

Typically these systems are studied in the adiabatic 
limit where the rate fl at which the driving force F = Qt 
is increased, is taken to zero allowing for a separation 
of the time scale of the fast local relaxation (given by 
the avalanche durations) from the scale set by the slow 
external driving (given by the rate fi). For experiments 
and commercial applications it is important to under- 
stand the effect of a finite but small driving rate fl > 
on the pulse statistics. In this letter we study this ques- 
tion, specifically focusing on the experimentally readily 
accessible BN in magnets as the primary example. Our 
results are fairly general in nature and expected to be ap- 



plicable to a much larger class of systems with crackling 
noise. 

A number of models have been used to study crack- 
ling noise ranging from interface depinning models and 
sandpile models, to models with many interacting inter- 
faces due to domain nucleation. While we believe the 
main results we report in this letter will apply to many 
approaches in the slow sweep rate regime we use a vari- 
ant of the zero temperature, nonequilibrium random field 
Ising model (ztneRFIM) to illustrate our results. The zt- 
neRFIM is a simple model for hysteresis and avalanches 
(BN) in magnets ||, ||. Neglecting the microscopic de- 
tails it assumes a hypercubic lattice of spins Si = ±1 with 
the Hamiltonian: 

where the first term represents the ferromagnetically 
coupled nearest neighbors (with coupling J > 0), 
H(t) = fit is the external applied field, and hi is 
the random magnetic field with Gaussian distribution 
p(hi) = cxp(-/i|/(2i? 2 ))/(i?v / 27r). hi models uncorre- 
cted quenched disorder. Other types of disorder, such 
as random anisotropies, and random ferromagnetic bonds 
are expected to be in the same universality class J9[ fl0| . 
We study the model at zero temperature, far from 
thermal equilibrium to model systems in which the ther- 
mal relaxation time scale is much larger than experimen- 
tal time scales. One sweep in our simulation begins at 
H = — oo causing all spins to be down. We then raise the 
external field slowly and flip each spin when its local 
effective field 

hf f = H(t) + J2^ + ^i (2) 

<*> 

changes from negative to positive. (Here (i) stands for 
the nearest neighbors of site i.) A single spin flip takes 
a microscopic time St = 1. It can trigger an avalanche 
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of spin flips propagated through the ferromagnetic inter- 
actions |(| . In the adiabatic limit O — > 0, the external 
magnetic field is kept constant during each avalanche, it 
is only increased in between avalanches, until the next 
avalanche is triggered. The avalanche time profile V(t), 
which is the number of spins flipped during the time in- 
terval of length St at time t is proportional to the voltage 
time profile V(t) measured in the induction coil in a cor- 
responding BN experiment during a BN pulse. 

For the adiabatic case it has been shown || that near 
a critical amount of disorder R — R c and a critical 
value of the magnetic field H = H c the avalanche size 
and duration distributions scale on long length scales 
as D S (S) ~ S- T V s {Sh 1 ^ a ^\h/r' 3S ) and D t (T) ~ 
T- a V t {Th uz /^ 5 \h/r fjS ), respectively, where r = \(R- 
R c )/R\ and h = \(H - H c )/H\, V s and V t are uni- 
versal scaling functions and a, v, z, /3, 5, a, r are all 
universal exponents that are known for the adiabatic 
case and are not all independent ||. The scaling rela- 
tion between the size and duration, T ~ S avz relates 
a to t through a — + 1. Numerical simulations 
render r = 1.60 ± 0.06 and a = 2.05 ± 0.05 for this 
model [jy]. Away from the critical point the correlation 
length £ ~ r _l/ S(/i/r' 3 ' 5 ) is finite, as is reflected by a cut- 
off to the power law scaling regime at a largest avalanche 
S* £ 1 / <T " and a largest duration T* ~ £ z . (Here S is 
another universal scaling function.) 

At finite O new avalanches can be triggered by the in- 
creasing external field before a previous one is completed. 
In this regime pulses observed in the time series V(t) do 
not necessarily correspond to the avalanches found in the 
adiabatic limit. (We stress here the subtle but necessary 
distinction between "avalanche" and "pulse".) Since there 
is no mechanism in the model for the sweep rate to effect 
the local dynamics of a spin flip (a spin flips in one time 
step if its local field goes from negative to positive regard- 
less of how large the positive local field is) the observed 
differences between pulses for finite fi, and avalanches 
for zero CI, are a result of the simultaneous propagation 
of avalanches as well as the possible internal changes to 
avalanches induced by spatial overlap of simultaneously 
propagating avalanches at finite ft. For comparison with 
experiments it is important to understand how the pulse 
size and duration distributions are affected by temporal 
and spatial overlap of the avalanches as ft is increased 
from zero. 

A number of experimental results have been reported 
in the literature regarding sweep rate effects on scaling 
exponents of pulse statistics [l^]. For BN in soft 
magnets without applied stress and with a(fl = 0) = 2 it 
has been reported that a and r exhibit a linear sweep rate 
dependence on fl as fl is increased from zero, whereas 
experiments with applied stress give a(Q = 0) < 2 and 
at finite Q the a and r exponents do not depend on fi. 
In the following we will show that these results can, in 
fact, be generalized and explained by straight forward 
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FIG. 1: Simulated BN voltage V (number of spins flipping 
per St) as a function of applied external magnetic field H ex t(t) 
(measured in units of coupling J), obtained using the ztneR- 
FIM. At the bottom one can observe separate pulses. As 
sweep rate is increased the pulses merge and become bigger 
until no separate pulses can be observed. 

temporal superposition of the adiabatic avalanches as ft 
is increased. 

There are two important threshold sweep rates, O t and 
fl s that mark the transitions to avalanches overlapping 
in time (fl > fit), and in space and time (ft > ft s ): 

(1) The lowest sweep rate Q = Oj at which on av- 
erage successive avalanches overlap in time so that one 
obtains a macroscopic uninterrupted pulse (see Fig [[]), is 
called the upper bound of "slow"(note that for f2 > Q t 
the number of pulses is greatly reduced making it diffi- 
cult to obtain pulse statistics at all). To estimate fl t we 
set the increase in H during the mean avalanche duration 
(T) equal to the mean separation ^p^y (in H) between 
successive avalanche seeds, and find 

a(H)(T) =i a(H)(T) ) l?' (3) 

where a(H) = a(H)L d is the non-universal number of 
avalanche seeds per unit field increase, and a(H) is inde- 
pendent of the system size L d (here L is the linear system 
size and d is the dimension of the system) . In general it 
is expected that fit will be proportional to l/n where n 
is the number of spins in principle available to nucleate 
an avalanche (volume in our model, and surface for sin- 
gle domain wall propagation). Eq.(||) implies that in the 
thermodynamic limit any finite sweep rate f2 results in 
the loss of distinct pulses, since all avalanches merge into 
only one big pulse. 

(2) The lowest sweep rate fl s at which successive 
avalanches tend to overlap not only in time, but also 
space can be estimated as follows: at fl s the expected 
number of new avalanche seeds triggered by the increas- 
ing external magnetic field within the volume of a large 
propagating avalanche of adiabatic size S, during its du- 
ration T, is a finite fraction of S. If r > 2 (i.e. (S) is not 
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divergent) £l s ~ ^_ ^> Q t (f or large enough L). If r < 2, 
in terms of the adiabatic quantities, we find 191 



Q g „ S £(t-2-<»/z)/(«0 _ ^ 

For large enough system sizes L, for finite correlation 
length £, one again obtains fl s ^> Therefore, the re- 
sulting time series V(t) in a BN experiment in the "slow" 
sweep rate regime < fit is simply a superposition of 
voltages produced by the individual simultaneously prop- 
agating avalanches. 

Superimposing power law distributed pulses yields re- 
sults consistent with the experiments. The reason for this 
requires discussion: on the one hand, the time between 
successive avalanche nucleations decreases for increasing 
sweep rate fi. On the other hand the duration of the 
avalanche results from the dynamics of the material and, 
in the slow sweep rate regime, does not depend on the fi. 
This is true for all avalanche models we have considered. 
The resulting pulse distributions are affected in two ways: 
(A) an avalanche of a certain size can appear to take 
longer, since a second avalanche starts before the first 
one ends, and runs longer ("swelling"). (B) Avalanches 
can completely disappear ("absorption") from the distri- 
butions due to their time interval of propagation lying 
entirely within that of another avalanche. 

(A) To illustrate the relationship between swelling 
and the adiabatic exponents we note that the largest 
avalanche nucleated during the propagation of another 
avalanche of duration T (with a > 1) is given by- 

rW a ~(a(fl)Tfi)^i (5) 

If a > 2 the avalanche of duration T will not swell ap- 
preciably since Tmax ~ T e where e < 1 and the new 

(2) 

"swollen" pulse is, at its largest, T + T ma x ~ T for large 

(2) 

T. For a < 2, however T m Lx > T is no longer negligible: 
by iteration the process runs away for very small sweep 
rate (for fi > fit ~ L d+}(2- a ) , see eqn(3) for a < 2). For 
a = 2 there is a linear increase in the duration that is 
T — > "fT . In this case, as we see below, absorption leads 
to a linear change of a with the sweep rate. 

The discussion for the pulse size distribution runs anal- 
ogously: The largest secondary avalanche triggered by 
the increasing field during a primary avalanche of dura- 
tion T has size (for r > 1) 

5il~( a (ff)Tfl)^ (6) 

The sweep rate significantly affects the distribution only 
if SmL > S. This occurs if > 1, i.e. if a < 2 in 
agreement with the above argument for the duration dis- 
tribution. If a < 2 the effects of sweep rate on the size 
distribution can be seen in |Lq] which we would argue 



results not in a decrease in the exponents, as suggested, 
but in a decreasing region S < S* over which the distri- 
bution can be described by a power law at all (see |15|| 

for more details) with S* ~ fi f^ 1 " 2 -") . 

(B) The relationship between absorption and the adi- 
abatic exponents is similar. As fi is increased from zero, 
note, that for an avalanche to maintain its identity (in 
terms of size and duration) it must not be absorbed 
into (i.e. propagate simultaneously with) larger ones. 
We assume the avalanches are uncorrelated in amplitude 
and have a nucleation density which can be described 
by a smoothly varying function of h and the probability 
T>(T, fi) that a certain avalanche of duration T remains 
a pulse of duration T as fi is increased is given by 

v(T,n) = D(T,n^o)- J] (l-pfT'.n-o)), (7) 

T'>T 

where P(T', fi — » 0) is the probability that the avalanche 
of duration T is absorbed in an avalanche of duration 
T', A detailed analysis of both swelling and absorption 
is given elsewhere . One finds that the effect of a slow 
but finite sweep rate depends on the adiabatic value of 
the exponent a: (1) For a(fi = 0) > 2, on long length and 
time scales, we expect no change in the scaling exponents 
due to the introduction of a slow sweep rate fi < fi t ■ (2) 
For a(fi = 0) = 2 we find a linear dependence of some of 
the exponents on the (slow) sweep rate: 

a(Cl) = 2~ K n r(Q) = - - — (8) 

where k is a non-universal constant. This is in agree- 
ment with the results found earlier for the ABBM model 
||l~7| , |T§| , |l9| , a single domain wall model , as well as for 
experiments in soft magnets with BN due to the propa- 
gation of parallel domain walls Jl3) : all of these systems 
have a(Cl = 0) = 2. (3) For a(fi = 0) < 2 the pulse 
perspective is lost even for very small sweep rates since 
temporal overlap takes over (leading to a large runaway 
pulse). For finite (but large) system size, in the slow 
sweep rate regime Q < fl t , so that one can actually still 
observe distinct pulses, these pulses will be distributed 
according to the adiabatic scaling forms, since in this 
case they tend to be very sparse on the time axis. This 
result also agrees with experimental results of Durin and 
Zapperi on soft magnets with stress induced anisotropics 
II- 

Numerically in the 3-d ztneRFIM we find that a(Q) 
appears to have a linear dependence on Cl since the zt- 
neRFIM has a = 2.05 ± .05, i.e. very close to 2 (see fig. 
||). We would expect, however, that if a>2, for very large 
durations (roughly T > 10" where a = 2 + e , for our 
problem e = .05) we would find, in the large system size 
limit, the finite but slow sweep rate to leave the exponent 
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batic limit. This conclusion is borne out by our numerical 
work on the ztneRFIM where we know that the pulses 
are simply a superposition of the adiabatic avalanches. 
We also provide criteria for determining for which sweep 
rates pulse statistics provide insight into the adiabatic 
collective behavior and when spectral methods will serve 
better. 
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PIG. 2: Distribution of avalanche durations T (number of 
time steps St = 1) at H ~ H c and R ~ R c for different sweep 
rates for the ztneRFIM in d — 3 dimensions. Inset shows 
linear sweep rate dependence of the exponent a taken from 
power law fits to the duration distributions. The simulated 
systems have L d = 500 3 spins. The data has been averaged 
over 20 disorder realizations. 



a unchanged. We estimate the system size required to 
bring the exponent to within the error bars (i.e. ±.05) of 
the adiabatic value for a = 2.05 and we find the system 
size must be orders of magnitude larger than the largest 
system we can simulate (500 3 ) (see for details of the 
base algorithm). It should be easier to observe the expo- 
nent independence in models with a significantly greater 
than 2. 

Conveniently it can be shown that the high fre- 
quency scaling behavior of the power spectra (PS) of the 
the BN voltage time series is independent of the sweep 
rate for Q, < fl s . This is a direct consequence of the time 
series being a mere superposition of the uncorrelated adi- 
abatic avalanches. For Poissonian distributed avalanche 
nucleation events (as in the ztneRFIM) the PS of the 
time series is a simple superposition of the PS of the in- 
dividual pulses ||22] , 23 1 and we note no difference in the 



scaling even in the lowest frequency of the PS. Since we 
are able to directly extract many of the exponents asso- 
ciated with the adiabatic avalanches (z,i/z.cti/2,§=£ 

\ ' ' 1 (XUZ 1 (7VZ ' 

from the PS regardless of sweep rate, collapses of the 
PS provide a useful tool for probing the collective be- 
havior without needing to address temporal overlap of 
avalanches. For O > f2 t where pulse statistics are lost 
and avalanche interactions become important as seen in 
|p4| one must use spectral tools. 

In summary we have found that the experimental 
observations of sweep rate dependence/independence of 
pulse duration and size distribution exponents can be 
understood as a consequence of temporal superposition 
of the avalanches that would be observed in the adia- 
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